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Abstract 

We address two finite size effects in perpendicular transport through multilayers of ferromag- 
netic (F) and normal metal (N) layers: (i) the transport properties depend on the magnetic layer 
thickness when of the order or thinner than the spin-flip diffusion length and (ii) magnetic layers 
with thickness approaching the magnetic coherence length become transparent for spin currents 
polarized perpendicular to the magnetization. We use magnetoelectronic circuit theory to investi- 
gate both effects on angular magnetoresistance (aMR) and spin transfer torque in perpendicular 
spin valves. We analyze recent aMR experiments to determine the spin-flip diffusion length in the 
ferromagnet (Py) as well as the Py|Co interface spin-mixing conductance and propose a method to 
measure the ferromagnetic coherence length. 

PACS numbers: 
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I. INTRODUCTION 



Since the discovery of the giant magnetoresistance (GMR)^ electron transport in mag- 
netic metallic heterostructures has been studied intensively and with considerable progress. 
The field developed from studies of large area multilayers of ferromagnetic (F) and normal 
metals (N) in which the current flows in the plane of the interfaces (CIP) to nanostructures 
with current perpendicular to the planes (CPP).- Current-induced magnetization excitation 
has been predicted for perpendicular F|N|F spin valves^ and subsequently observed. 
In these experiments applied currents excite a spin accumulation in the normal metal spacer 
that exerts a torque on the ferromagnets. When this torque overcomes the damping, the 
magnetization starts to precess coherently, possibly leading to a complete magnetization 
reversal.- By fits of the parameters of the diffusion equation^ to a wealth of experimental 
data of the GMR in CPP structure, the spin-dependent interface and bulk material re- 
sistances of the most important transition metal combinations are well known by now.- 1 ^ 
First-principles calculations in general agree well with the experimental values.— Also in 
view of possible applications for switching purposes in magnetic random access memories, 
a comparably accurate modeling of the spin torque as a function of material combinations 
and applied bias is desirable. 

Physically, the spin-transfer torque is a consequence of angular momentum conservation 
when a spin current polarized transverse to the magnetization direction is absorbed at the 
magnetic interface.— The transverse spin current can penetrate the ferromagnet up to a 
skin depth equal to the ferromagnetic coherence length A c = n/ | — kf\ . In transition 
metals A c is much smaller than all other length scales such as spin-diffusion length or mean- 
free path.— dtiJJi When the ferromagnetic layer thickness dp ^> A c the spin-transfer torque 
is a pure interface property governed by the so-called spin-mixing conductance,— which is 
accessible to first principles calculations.— 

An excellent method to measure the torque and mixing conductance is the normalized 
angular magnetoresistance (aMR) of perpendicular F|N|F spin valves^k^ 

R(9)-R(0) 

where R (9) is electric resistance when the two magnetizations are rotated by an angle 9 
with respect to each other. Deviations of the aMR as a function of cos 9 from a straight 
line are proof of a finite mixing conductance.— Systematic new measurements of the aMR 



have been carried out recently by Urazhdin et alM- on Permalloy(Py)|Cu spin valves as a 
function of the Py thicknesses. 

Interesting effects such as non-monotonic aMR, change of sign of the spin-transfer torque 
and strongly reduced critical currents for magnetization reversal have been predicted for 
asymmetric spin valves.— Asymmetry here means that the two ferromagnets in the spin 
valve are not equivalent for spin transport. This can be achieved by different thicknesses 
of the magnetically active regions of otherwise identical ferromagnetic contacts, but only 
when the spin-flip diffusion length in the ferromagnet lf d is of the order or larger than one 
of the magnetic layer thicknesses. The magnetically soft Py is the material of choice, but its 
spin-flip diffusion length is only lf d = 5 nm.— Urazhdin et a/.— investigated spin valves with 
ultrathin cIf < ifdi which means that the analysis of these experiments requires solution of 
the spin and charge diffusion equation in the ferromagnet. 

Detailed calculations for transition metals^ 1 ^ confirm that a transverse spin current can 
penetrate the ferromagnet over distances < 1 nm as a consequence of incomplete destruc- 
tive quantum interference. Urazhdin et al. investigated spin valves with Py layers of such 
thicknesses, claiming to observe an effect of this transverse component on the aMR. In weak 
ferromagnets like CuNi or PdNi alloys in which A c may become larger than the scattering 
mean-free path, the transverse component of spin current and accumulation can be treated 
semiclassically.— It is shown below that an effective conductance parameter ("mixing trans- 
mission") can be introduced to parametrize transport in both regimes. 

In this paper we treat the size effects related to dp < lf d (Section II) and dp % A c 
(Section III) (but A c much smaller than the spin diffusion length). In Section II we apply 
magnetoelectronic circuit theory^ combined with the diffusion equation to the F|N|F|N 
spin valves studied by Urazhdin et al.. We demonstrate that the angular magnetoresistance 
provides a direct measure for the mixing conductance^ and find that the non-monotonicity 
in the aMR is indeed caused by the asymmetry as predicted. For F|N|F|N|F structures, that 
are also of interest because of their increased spin torque,—^ we obtain several analytical 
results. The approach from Section II is generalized in Section III allowing us to treat 
ultrathin ferromagnetic layers or weak ferromagnets.—^ 1 ^ We find that there should be no 
measurable effects of A c on the aMR in F|N|F|N structures, but predict that the torque 
acting on the thin layer is modified. We proceed to conclude that the coherence length 
should be observable in the aMR of F|N|F|N|F structures. Finally, we propose a set-up to 



3 



measure the ferromagnetic coherence length in a three-terminal device. 



II. MAGNETOELECTRONIC CIRCUIT THEORY AND DIFFUSION EQUATION 
FOR SPIN VALVES 

In this Section we assume that A c <C dp. In Part A we recapitulate some old results: 
the magnetoelectronic circuit theory for spin valves, with emphasis on the inclusion of the 
spin-flip diffusion in the ferromagnetic layers when the ferromagnetic layer thickness dp is 
of the same order as the spin-flip diffusion length in the ferromagnet lf d . In Part B we apply 
these results to recent experiments by Urazhdin et al. in which we can disregard spin-flip in 
the Cu spacers. In Part C we present new results for symmetric F|N|F|N|F structures. 



A. Magnetoelectronic circuit theory and diffusion equation 

Magnetoelectronic circuit theory^ has been designed to describe charge and spin trans- 
port in disordered or chaotic multi-terminal ferromagnet-normal metal hybrid systems with 
non-collinear magnetizations. The material parameters of the theory are the bulk and inter- 
face spin-dependent conductances, as well a the so-called interface spin-mixing conductance 
Gfj_. For spin valves, circuit theory can be shown to be equivalent to a diffusion equation 
when ImG^i ~ 0, which is usually the case for intermetallic interfaces.— When the thickness 
of the ferromagnetic metal layer d 3> lf d , the layer bulk resistance can be effectively replaced 
by that of a magnetically active region close to the interface of thickness lf d . When connected 
to a reservoir or other type of spin sink, the effective thickness becomes lf d ta.nh(dp/lf d ).— 

The aMR for general N|F|N|F|N structures with ImG|| = as derived previously 2 ^ reads 

Rn(Ri- +aR 2 -) 2 + (1 -a 2 )(R 2 R 2 + R 2 (R 1 + R n )) 
m = R n + R1 + R2- ( " } - H - 2 2 ~ { (2) 

with a = cos 9, 4Ri (2 ) = l/Gi^n + l / G K?)l ~ 2R n> 4 -Ri(2)- = 1/^1(2)1 ~ 1 / G i(2)|, Pi(2) = 
Ri(2)- / R\(2), and 2R^ = + l/G^Ti) where Gi(2)t an d ^?i(2)j, are conductances of the 

left (right) ferromagnet including the left (right) normal layer, Gi^ and G%\\, are mixing 
conductances of the middle normal metal with adjacent ferromagnet interfaces as shown in 
Fig. 1. The torques felt by first and second ferromagnetic layer become 

h 1 + RtdRx - aPJP 2 

Tl / J ° = o~ TTjTr /B Vi -a-g m\ 2 m i x m2 x mi 3 

2e (1 + R n /Ri)[l + R u /R2) - a 2 



Figure 1: Definition of conductances Gi(2)m) an d mixing conductances G 1 ( 2 )fj. f° r N|F|N|F|N 
structure. 
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2e(l + R n /R 1 )(l + R n /R 2 )-a 
When we approximate the mixing conductance 1/Rfi by the Sharvin conductance of the 
normal metal, Eqs. (}3f4l) coincide with the expressions in Ref. H for asymmetric N|F|N|F|N 



spin valves with A 



L(R) = 2R m/ R U + !> P L(R) A L(R) 



2Ri(2)-/R 



B. Extraction of the mixing conductance from experiment and asymmetric spin 
valves 

Most material parameters in circuit theory are those of the two-channel resistor model. 
They can be determined for the collinear magnetic configurations, i.e. via the (binary) GMR. 
The only additional parameters for the non-collinear transport are the interface mixing 
conductances CrL, , assumed here to be real. These can be found from a single parameter fit 
of the experimental aMR or from band structure calculations. A symmetric F|N|F structure 
is most suitable to carry out this program. The thus obtained G?*, should be transferable 
to other (asymmetric) structures grown by equivalent techniques. Urazhdin et al. fitted 
their experimental results for the normalized aMR by the simple formula^ that follows from 
circuit theory^ 

aMR(g)= /1 1 ~ C °'f , (5) 
v ; X(l + cos#) + 2 



For symmetric junctions we identify x = 2R/R^i (see Eq. (18) in Ref. 124 ). 

Urazhdin et al& used the structures Nb(150)Cu(20)FeMn(8)Py(rfi)Cu(10)Py(rf 2 )Cu(20)Nb(150), 
where the numbers in brackets are the thicknesses in nm. The exchange bias antiferromag- 
net FeMn is treated as a perfect spin sink, which means that the effective thickness of the 
left Py layer becomes lf d t&nh(di/lf d ) = 0.8l s d (d\ = 6 nm, lf d = 5.5 nm). Note that this 
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device is not exactly symmetric when d 2 3> lf d as di is not much larger than lf d . but the 
calculated deviations from the fitted mixing resistances are smaller than the experimental 
error bars. When we replace d\ by lf d and d 2 3> lf d the sample is symmetric and the aMR 
is well represented by Eq. (03) with \ = 1.96 (see Fig. 2). 23 

We can use the measured value of x to derive the mixing conductance l/(ARn) of an 
interface with area A by = 2Rj\. For comparison with first principles calculations for 
point contacts based on the scattering theory of transport,— the Sharvin resistance of the 
normal metal should be added^ 2 - 

R p {{ = R-u + R s h 

Using the notation: 

AR = l S dp*p y + AR* PyCu - AR n /2 
AR_ = lsdP*Pyf3 P y + AR* PyCu ^ 

we may substitute the well established material parameters for bulk and interface resistances 
of Cu|PyM l sdP * Py = 1.4 fftm 2 , l sd = 5.5 nm, AR* PyCu = 0.5 ffim 2 , f3 Py = 0.7, 7 = 0.7, 
disregarding the small bulk resistance of Cu which led us to AR^ = 1.3 fflm 2 and t] = 1.49. 
This value of the mixing resistance is larger than the Sharvin resistance AR s h — 1/G — 
0.878 fflm 2 of Cu used by Xiao et al— The point-contact mixing resistance of the Cu|Py 
interface that should be compared with band structure calculations is AR*?' = 2.2 ff2m 2 , 



somewhat smaller than that found in Ref. |22| (2.56 ffhn ). Both results are close to the 
band structure calculations 1 - of the point-contact mixing resistance for the disordered Cu|Co 
interface (2.4 ff2m 2 ). 

In Fig. 2 we compare plots of Eq. (EJ) with experimental aMR curves for symmetric and 
asymmetric F|N|F|N|S multilayers,— identifying the following relations between parameters: 

ARi = l sd ppy + AR* PyCu - AR n /2 

AR 2 = dip* Py + AR* PyCu + ARpym — AR-^/2 

AR X _ = l sd p* Py P Py + AR PyCu ^ 

AR 2 _ = d 2 p* Py [3 Py + 2AR* PyCul . 

We assume that the spin current into the superconductor vanishes. The resistance between 
the right ferromagnet and the right reservoir was taken to be ARp y Mb = 5 ff2m 2 . This is 
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Figure 2: aMR of the F|N|F|N structure for two thicknesses of the right ferromagnetic layer d<i = 
0.27l sc i, 2.2l s d (d\ > l s d). The filled (large d?) and open (small d?) squares are the experimental 
data. 23 
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Figure 3: Angular magnetoresistance — K(0) of the F|N|F|N structure for different thicknesses 
of the right ferromagnetic layer d% = 0.27/ S( 2, 0.5/ s rf, 2l S( i, 2.5l s d and oo (starting from the lower 
curve respectively (di 3> lf d , ARp y Ni, = 3 ffim 2 )). 
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larger than the AR PyNb = 3 ff2m 2 reported in Ref. |34|, but gives better agreement with the 
experiment. We observe good fits in Fig. 2 nicely reproducing the non-monotonic behavior 
around zero angle. 

In Fig. 3 we plot the angular magnetoresistance for different thicknesses of the right 
Py layer, all relative to the parallel configuration, but not normalized to a relative scale 
as above. The lower curve was obtained from Eq. ([2]), the others were calculated numeri- 
cally solving the bulk layer spin-diffusion equation in the ferromagnet. The non-monotonic 
angular magnetoresistance disappears when the right ferromagnetic layer becomes thicker 
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Figure 4: Angular magnetoresistance — K(0) and spin-transfer torque on the left ferromagnet 
for the F|N|F|N structure with right F-layer thickness di = 0.27/ s rf {ARp y Nb — ► oo, tq = I^hjle). 




and therefore the sample more symmetric. For the set of parameters in Fig. 3 the non- 
monotonic behavior is rather weak but with circuit theory we can readily propose samples 
that maximize the effect. The minimum of the angular magnetoresistance Eq. j2|) at finite 
9 1 that coincides with a zero of the spin-transfer torque on the left ferromagnet^ 1 ^ 

(R n + Ri)R 2 - 



COs6>! 



(6) 



R1-R2 

To observe the effect clearly, cos^i should be small, which can be achieved by increasing -R2+, 
e.g. by the resistance of the right-most normal metal (within the spin-flip diffusion length). 
In Fig. 4 we plot the angular magnetoresistance Eq. ([2j) and the spin-transfer torque on the 
left ferromagnet Eqs. ([3]) when the resistance of the right contact is AR PyNb — > 00. 



C. Analysis of symmetric F|N|F|N|F structures 

Our approach offers analytic results for symmetric F|N|F|N|F structures when the outer 
layers are thicker than lf d . In Fig. 5 we plot the angular magnetoresistance when the 
magnetizations of the outer layers are kept parallel for material parameters that are the 
same as above and close to set-up B from Ref. 23]. When the middle layer thickness d 3 ^> lf d 
the angular magnetoresistance is equal to that of two symmetric F|N|F structures in series. 
The analytical formula for the angular magnetoresistance in the regime d 3 -C lf d is presented 
in Appendix A. For d 3 > 0.3/^ we cannot disregard spin flip in the middle layer and compute 
the resistances numerically. 

A symmetric F|N|F|N|F setup with antiparallel outer layers can increase the torque.— 
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Figure 5: Angular magnetoresistance 9?(#) — 9?(0) of the F|N|F|N|F structure for the middle F layer 
thicknesses d = 0.27l s d, 0-5l s d, ^hd, 2.5/ s d and oo (starting from the lower curve, respectively). The 
parallel resistance is subtracted. 




Enhancement by a factor of 2 was reported by S. Nakamura et al.—. This result can be ob- 
tained from the magnetoelectronic circuit theory (see also S. Nakamura et al., unpublished). 
With a current bias Jo, assuming d 3 <C lf d , we derived a simple formula (note the similarity 
with the torque on the base contact of the three-terminal spin flip transistor—): 

h 2 J R_|sin#| 

T// ° " TeR n + R^6' (7) 
without invoking the parameters of the middle layer. When d 3 ^> l S( i we can divide system 
into two F|N|F spin valves in series. Taking into account Eq. (3), the torque can be written 
down immediately: 

t/Iq = T FNF {6)/I + T FNF (<K - 9)/I (8) 

H -R_ | sin 9 1 h i?_ | sin 9 \ 

~ 2eR u + R{l + cos 9) + Ye R n + R(l - cos 9) 

In Fig. 6 we plot results of these two analytic formulas as well as results of numeric calcula- 
tions for the case d = 0.8Z s d- Note that these curves are symmetric with respect to 9 = tt/2. 
By the dashed line we plot the torque for the corresponding symmetric F|N|F structure. 



III. COHERENT REGIME 

Urazhdin et a/.— 's intentions to search for coherence effects in ultrathin magnetic layers 
encouraged us to study the regime d F < A c . In this Section we formulate the magnetoelec- 
tronic circuit theory that includes coherence effects in this regime in two and three terminal 
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Figure 6: The spin-transfer torque on the middle ferromagnet for the F|N|F|N|F structure for the 
thickness of the middle layer d = 0.27Z st j, 0.8l s d and Wl s d starting from the lower curve respectively 
(by bold line), the same for the corresponding symmetric F|N|F structure (by dashed line), tq = 
I h/2e. 




multilayer structures. Since A c is only a couple of monolayers, we are allowed to disregard 
spin-flip and diffuse scattering in the ferromagnetic material bulk layers. 

A. Extended magnetoelectronic circuit theory 

We consider an N1|F|N2 circuit element, choosing the normal metals as nodes with a 
possibly non-collinear spin accumulation and the entire F layer including the interfaces 
as resistive element (see Fig. 7). This allows us to treat the ferromagnet fully quantum 
mechanically by scattering theory. The current through the ferromagnet depends on the 
potential drop between and the spin accumulation in each of the normal metal nodes. Spin 
I s and charge Jo currents can conveniently expressed as 2 x 2-matrices in Pauli spin space 
/ = (l/o + <7 ■ Is)/2, where a is the vector or Pauli spin matrices and 1 the 2x2 unit matrix. 
On the normal metal side^ in the region 2 

f= e^^nm^ny _ ^ 2 + ^ 2(pm)t] } (g) 
nm 

where r mn is the spin dependent reflection coefficient for electrons reflected from channel 

^-mn 

n into channel m in the node 2, t' is the spin dependent transmission coefficient for 
electrons transmitted from channel n in the node 1 into channel m in the node 2 and S nm is 
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Figure 7: A contact through a thin ferromagnet between two normal metals nodes. The current is 
evaluated in the node 2. 
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the Kronecker delta symbol. 

^~-mn 

In the absence of spin nip processes the matrices r mn and r should be diagonal in 
spin space provided the axis z is parallel to the magnetization of the ferromagnet (we are 
free to chose this frame reference as it is more convenient). Expressing the spin-dependent 
distribution matrices in nodes 1 and 2 via Pauli matrices; f N = l/f + crff and the unit 
vector m 2 parallel to the axis z we obtain for spin and charge currents in the node N2: 

Jo = (G T + GJAf" + (G T - GJAff • m 2 , (10) 



m 2 [(G T - GJA/f + (G T + GJAff ] 



-2(m, x ff 2 x m z ) ReGjf 2 |F + 2(ff 2 x m z ) ImG^f 2 ' 7 

4f |JV2 -2(ff 1 xm z )ImGf i 



+2(m z x ff 1 x m z ) ReGlf m - 2(ff 1 x mj ImG^ 2 . (11) 



where A/f = if 1 - and Aff = ff 1 - ff 2 . This agrees with the result of Ref. 
except for the terms involving the mixing transmission^ 1 *^ 2 - 
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2 2 
n rN2\F _ e V^/r nm ( nm\*\ n tNl\N2 _ & ST~^ . Inm ( ,lnm\* 

nm nm 

The torque acting on the magnetization through the interface adjacent to N2 is the transverse 
component of the spin current flowing into the ferromagnet: 

r 2 = -2(m 2 x ff 2 x m t ) ReG^f 2 |F + 2(ff 2 x m z ) ImGjf 2|F 

+2(m 2 x ff 1 x m 2 ) ReG^ 2 - 2(ff 1 x m 2 ) ImG^ 2 . (12) 



When two opposite direction of the magnetization M and — M are equivalent for the 

tN: 

n 



transport, we obtain Glf 2 = G 1 ^ 1 as a consequence of time reversibility. This condition 
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should hold in most cases (e.g. Stoner model is isotropic in spin space). The mixing 
transmission describes the part of the transverse spin current that is not absorbed by the 
ferromagnet and vanishes when the ferromagnetic layer is thicker than the ferromagnetic 
coherence length A c .— It is complex, its modulus representing the transmission probability 
and the phase of the rotation of the transverse spin current by the ferromagnetic exchange 
field. First-principles calculations of G*. and Gt, have been carried out by Zwierzycki et 
alrZ- showing small variation of the first and non-vanishing value of the second when the 
ferromagnetic layer becomes of the order of several monolayers. 



B. Observation of ferromagnetic coherence in transport experiments 

In this section we address coherence effects due to the transmission of transverse spin 
currents through ultrathin ferromagnetic layers or weak ferromagnets. These effects should 
be observable in Py structures when dp < 1.5 nm. Band structure calculations show that in 
Cu|Co|Cu structures the mixing transmission can easily reach Gtt ~ 0.1 (Gj + GJ for such 
thicknesses.— 

We may draw an important conclusion from the extended magnetoelectronic circuit the- 
ory applied to general (asymmetric) N1|F1|N2|F2|N3 structures: when the nodes are chosen 
in the middle normal metal and in the outer normal metals at the points that connect to 
the baths, a possibly finite mixing transmission completely drops out of the charge trans- 
port equations, i.e. the expressions remain exactly the same as those derived above for the 
N1|F|N2 structure. For example, the charge and spin currents from Nl (and similarly from 
N3) into N2 read 

I = (G T + G l )Afi Q + (G T - G l )A/j, s • m z (13) 



m. 



{G m\m _ G ™m )AfJ N + (G N1|N2 + G NHN2 )A/4 7 



-2(m z x M f x m z ) ReG^ 2|F1 + 2( M f x m z ) ImG^ 2|F1 

+2(m z x 1 x m z ) ReG^ 2 - 2(/xf 1 x m z ) ImG^ 2 , (14) 

where A/j,q = fi^ 1 —/!^ 2 and A/j,^ = n^ 1 —^ 2 describe the potential and spin accumulation 
drops between the left and the middle nodes. By conservation of spin and charge currents 
in the center node, expression for aMR can be derived. However, the mixing transmission 
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does not appear in Eqs. (1 1311 141) since there is no spin accumulation in the outer nodes 
(reservoirs). Ferromagnets thin enough to allow transmission of a transverse spin current 
can therefore not be distinguished from conventional ones in the aMR. Our conclusions 
therefore disagree with the claims of ferromagnetic coherence effects in aMR experiments on 
N|F|N|F|N structures by Urazhdin et al..— 

On the other hand, the torque on the thin ferromagnet F2 does change: 

t 2 = -2(m z x [is x m z ) Re(G n 1 - G n ) + 2(ti, x m z ) Im(G n 1 - G u ). 

A parameterization of the torque via a combination G^ 2 ' F2 — Gt 1 ^ 2 ^ 3 was found in Ref. 



13l by random matrix theory, which is equivalent with circuit theory when the number of 
transverse channels is large.^ However these authors did not discuss their results in the limit 
of thin ferromagnetic layers. When ImG'| ~ and ImG*, ~ 0, the torque r coh acting on 
the thin layer is modified from the incoherent expression r as: 

r«* = T(cr n -G t n )/cr n . (15) 

Naively one may expect that the reduced absorption of the transverse spin accumulation 
diminishes the torque, but this is not necessarily so (see Fig. 8). Since the mixing trans- 
mission may be negative, Eq. ffl~5l) shows that increased torques are possible. This can be 
understood as follows. A spin entering a ferromagnet will precess around an exchange field 
normal to its quantization axis. A negative mixing transmission ReGj^ < adds a phase 
factor corresponding to a rotation over an angle tt during transmission. The outgoing spin 
then has a polarization opposite to the incoming one. The magnetization torque, i.e. the 
difference between in and outgoing spin currents, consequently increases compared to the 
situation in which the incoming transverse spin is absorbed as in thick ferromagnetic layers. 

In contrast to N|F|N|F|N structures, we find that it is possible to observe Gl> in the aMR 
of F|N|F|N|F devices. We study here the dependence of the aMR on the mixing transmission 
in a Py based multilayer. In Fig. 9 we present the aMR for different mixing transmissions 
in the middle layer of thickness oIf = 0.27Z S ^- Unfortunately, it seems difficult to obtain 
quantitative values for the mixing transmission from experiments since the dependence of 
the aMR on Giy is rather weak. 

When the coherence length becomes larger than the scattering mean-free path, which can 
occur in weak ferromagnets like PdNi or CuNi, the transverse spin accumulation should be 
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Figure 8: The torque on the thin right layer of thickness d = 0.27l s d for F|N|F structure. The 
left layer has thickness d ^> lf d . The curves starting from the lower one respectively, Re(l/G*^) = 
5ffhn 2 , Im(l/G* n ) = oo; Re(l/G* n ) = oo, Im(l/G' n J = 5mm 2 ; Re(l/G* n ) = oo, Im(l/G* u ) = oo; 
Re(l/G* ) = -5mm 2 , Im(l/Gt ) = 00 (r = I h/2e). 
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Figure 9: aMR of F|N|F|N|F structure for the thickness of the middle layer d = 0.27/ S( ^. Outer layers 
are antiparallel with d ^> lf d . The curves starting from the lower one respectively, Re(l/Gi,) = 
-5mm 2 , Im(l/G* n ) = 00; Re(l/G* n ) = 00, Im(l/G^) = 5m m 2 ; Re(l/G* n ) = 00, Im(l/G* u ) = 
00; Re(l/Gl) = MQm 2 , Im(l/Gl) = 00. 
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treated by a diffusion equation.— The result can be parametrized again in terms of a mixing 
transmission, which can subsequently be used in our circuit theory. 

C. Three terminal device for observation of coherence effects 

Finally, we propose an experiment that should be more sensitive to ferromagnetic co- 
herence. We suggest the setup shown in Fig. 10 that is analogous to the spin-torque 
transistor^ and the magnetoelectronic spin-echo^ 2 - concepts. A current through the antipar- 
allel ferromagnets Fl and F2 excites a spin accumulation in the normal metal Nl. This spin 
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accumulation can transmit F3 only when its thickness is less than A c . In that case a spin 
accumulation is induced in the upper normal metal N2 that can be detected as a voltage 
depending on the magnetization angle 9 of the analyzing ferromagnet F4. We assume here 
that Nl is smaller than its spin-flip diffusion length (Cu is a good candidate with spin- 
diffusion lengths of up to a micron) such that the spin accumulation is constant under the 
contact to F3. Otherwise the signal at the ferromagnet F4 is diminished since part of the 
spin accumulation in Nl is lost due to spin-flip processes. 

When the Gi, of F3 is smaller than its and ferromagnet F4 is not too leaky for the 
spin current (e.g. connected via a tunnel junction) the spin accumulation in N2 can be 
found from Eqs. (1101)111) in terms of the spin accumulation in Nl: 

/ o N 



(ReG r n ) 2 + (lmG r n ) 2 
( 



Re G r n Re G t n + Im G r n Im G t n 



ReG r u lmG t n 



(16) 



ImG^ReG^y 



l/%ll 

ReG r u 





ReG 
ImG 
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where Eq. (|T7j) holds to a good approximation when the layer F3 is metallic. The spin 
accumulation is indeed coherently rotated by the exchange field in F3. The angle dependence 
of the potential in F4 is U{6) ~ fi sl P\Gi^\ cosO/ (ReGlA with maximum along fi S2l where 
P is the polarization of the contact N2|F4. 

When the G*^ of F3 is not smaller than Gl, (or the spin current leak into F4 is significant), 
the spin accumulation n si is affected by fi S2 and the final expressions are more complicated. 

An angle dependence of U (9) provides a direct proof of a finite mixing transmission. The 
ferromagnetic coherence length can be determined by repeating experiments for a number 
of layer thicknesses of F3. Such a direct experimental evidence should help to get a grip on 
this important parameter A c . 



IV. CONCLUSION 

In this paper we extracted the spin-mixing conductance of a Py|Cu interface from the ex- 
perimental data of Urazhdin et al. using material parameters measured independently by the 
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Figure 10: An experimental setup to observe the mixing transmission and measure the ferromagnetic 
coherence length in ferromagnet F3. The spin accumulation f2 S2 in the normal metal N2 is measured 
via the angular dependence of the potential U(9) of the ferromagnet F4 that is weakly coupled to 
N2. 
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MSU collaboration. We find good agreement with experiments on asymmetric F|N|F|N mul- 
tilayers, reproducing quantitatively the non-monotonic aMR that we predicted earlier.— 
Magnetoelectronic circuit theory together with the diffusion equation is a convenient tool 
for the data analysis when the spin-flip diffusion length in the ferromagnet is of the same 
order as the layer thickness. We suggest carrying out current-induced magnetization reversal 
experiments on samples that display the non-monotonic aMR since we predict anomalous 
magnetization trajectories due to a vanishing torque at finite magnetization angle.— We 
also study the effects of the finite ferromagnetic coherence length in ultrathin ferromagnetic 
films or weak ferromagnets. For this purpose the magnetoelectronic circuit theory is ex- 
tended to treat phase coherent transport in the ferromagnet. A coherence length that is 
larger than the ferromagnetic layer thickness does not modify the aMR of N|F|N|F|N struc- 
tures, but a small effect should exist in F|N|F|N|F structures. In contrast, the spin-transfer 
torque is affected more strongly and may even be increased by the spin-coherence when 
the exchange field rotates the transverse spin current polarization by the angle n. Finally, 
we propose a three-terminal device that should allow experimental determination of the 
ferromagnetic coherence length. 
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Appendix A: ANALYTICAL RESULTS FOR F|N|F|N|F STRUCTURE. 

The aMR of a F(T)|N|F(0)|N]F(T / j) CPP pillar can be described analytically when the 
thick outer layers are parallel or antiparallel, respectively. 



R{9) = 2(R n + R) + R M 

Rn( R li~ + 4g-(g- + R M -a)) + (2R M R 2 _ + R 2 M ^R)(i - a 2 ) 
(R n + R)(2R n + R M ) - RR M a 2 : 

2R 2 Jl-a 2 ) 



(Al) 



R{9) = 2{R n + R) + R M - 



R u + R(l-a 2 ) 
2R 2 M _(R n + R(l - a 2 )) 



(A2) 



(R n + R)(2R U + R M ) - RRm® 2 ' 

where a = cos 9, 4R + 2R^ = + ^ , 4_R_ = ^ — ^- for the outer layers. The mixing 
resistance for two interfaces adjacent to any normal metal Rn = (we assume all interfaces 
identical). For the middle layer 4Rm = ^ + — , ^Rm- = ^ — ^j- 
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